A = imread('watermarked.tif');

B = imread('cameraman.tif');

alpha1 = 0.98;

T = dFRT(256,alpha1);
F1 = T*double(A)*T;

% [mag_array loc] = sort(abs(F1(:)),'descend');
% 
% percent = .17;
% 
% start_pt = 1+floor(percent * length(loc));
% 
% extract1 = F1(loc(start_pt:start_pt))

F2 = T*double(B)*T;

F = F1-F2;

figure,surf(abs(F));

alpha2 = 0.93;

T = dFRT(256,alpha2);

X = T*F*T;

figure,imshow(mat2gray(abs(X)));

